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Modeling real-world systems plays a piv- 
otal role in their analysis and contributes 
to a better understanding of their behav- 
ior and performance. Classification, opti- 
mization, control, and pattern recognition 
problems rely heavily on modeling tech- 
niques. Such models can be categorized 
into three classes: white-box, black-box, 
and gray-box (Nelles, 2001). White-box 
models are fully derived from first prin- 
ciples, i.e., physical, chemical, biologi- 
cal, economical, etc. laws. All equations 
and parameters are determined from the- 
ory. Black-box models are based solely 
on experimental data, and their struc- 
ture and parameters are determined by 
experimental modeling. Building black- 
box models requires little or no prior 
knowledge of the system. Gray-box mod- 
els represent a compromise or combina- 
tion of white-box and black-box models 
(Nelles, 2001). 

In the modeling of highly nonlinear 
and complex phenomena, we may not 
have a good understanding of the pro- 
cesses, and thus black-box models may 
be our best (or even our only) choice. 
Artificial neural networks (ANNs) are one 
of the most powerful and popular tools 
for black-box modeling and are designed 
and inspired by real biological neural net- 
works. 

There has been an increasing interest 
in analyzing neurophysiology from a non- 
linear and chaotic systems viewpoint in 
recent years (Christini and Collins, 1995; 
Sarbadhikari and Chakrabarty, 2001; Korn 
and Faure, 2003; Hadaeghi et al., 2013; 
Jafari et al., 2013; Mattei, 2013). For 
example, although the famous Hodgkin 
and Huxley model (Hodgkin and Huxley, 



1952) has been the basis of almost all 
of the proposed models for neural firing, 
the Rose-Hindmarsh model (Hindmarsh 
and Rose, 1984) is known to be a more 
refined model since it can show different 
firing patterns, especially chaotic bursts of 
action potential, which enable a proper 
matching between this model behavior 
and experimental data. Another example 
of the observation of chaotic behavior in 
the nervous system is the period-doubling 
route to chaos in flicker vision (Crevier 
and Meister, 1998), which is the focus of 
this letter. 

Stimulation with periodic flashes of 
light is useful for distinguishing some 
disorders of the human visual system 
(Crevier and Meister, 1998). It has been 
shown by Crevier and Meister (1998) that 
during electroretinogram (ERG) record- 
ings of the visual system, period-doubling 
can occur. It is well-known that period- 
doubling occurs in nonlinear dynami- 
cal systems, and it is often associated 
with the onset of chaos. In one study 
(Crevier and Meister, 1998) the retina of 
a salamander was stimulated with a peri- 
odic square-wave flashes, and the ERG 
was recorded. The flash frequency was 
changed between zero and 30 Hz, while 
the contrast was constant. In another 
record, the contrast was changed while 
the frequency was fixed at 16 Hz. All the 
ERG signals were filtered at 1-1000 Hz. 
Using a common approach to obtain 
a discrete time series from a continu- 
ous recorded signal, successive local max- 
ima of the signal were extracted as a 
time series (Figure 1A). As shown in 
Figures 1B,C, both the parameters (flash 
frequency and contrast) have a great effect 



on the recorded ERG signals and cause 
bifurcations resulting in a period-doubling 
route to chaos. 

However, it is difficult to understand 
the exact relations between the parame- 
ters and their effects. In other words, it 
is not easy to build a white-box model 
that can regenerate the signals and dia- 
grams accurately. That may be because of 
the highly complex and nonlinear dynam- 
ics involved. We have used the ability of an 
ANN in learning highly nonlinear dynam- 
ics as a black-box model of this system. We 
used a four hidden layer feed-forward neu- 
ral network with (7/4/8/5) neurons in the 
layers (Figure ID) and nonlinear trans- 
fer functions hyperbolic tangent function 
that help the network learn the com- 
plex relationships between input and out- 
put. The activation function of the last 
layer of the network is linear transfer 
function. We used two parameters (con- 
trast and frequency) and three time delays 
{Xn-i, x n -2-> and Xn-i) as the inputs of 
the ANN to fit each data point of the 
time series (x n ) as the output of the 
network. 

As shown in Figures 1E,F, this model 
can generate bifurcation diagrams simi- 
lar to those obtained from real data. As 
the result, we believe that ANNs are pow- 
erful tools for modeling highly nonlin- 
ear behavior in the nervous system. We 
plan to construct ANN models in future 
work including extension to more cases 
and details, extension of the ideas in 
Hadaeghi et al. (2013) to patients with 
bipolar disorder, and extension of the ideas 
in Jafari et al. (2013) to patients with 
attention deficit hyperactivity disorder 
(ADHD). 
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FIGURE 1 | (A) One example of the local maxima of the ERG signals. (B) Real 
bifurcation diagram resulted from varying flash frequency. (C) Real bifurcation 
diagram resulted from varying contrast. (D) The structure of the ANNs that 



were used. (E) Artificial bifurcation diagram resulted from varying the flash 
frequency input in the ANN. (F) Artificial bifurcation diagram resulted from 
varying the contrast input in the ANN. 
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